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SITE ACCESSIBILITY AND CHARACTERISTIC VELOCITY REQUIREMENTS 
FOR DIRECT-DESCENT LUNAR LANDINGS 
by Vernon J. Weyers 
Lewis Research Center 

SUMMARY 

Direct -descent lunar landings are shown to be possible for landing sites located at 
any point on the lunar surface. The direct descent characteristic velocity requirement 
is shown as a function of landing site location for various combinations of descent engine 
specific impulse, translunar trip time, and initial thrust to Earth-weight ratio values. 

The landing site location is referenced to the normal impact point of the incoming 
trajectory. A means of estimating the normal impact point is provided for cases in which 
it is not known. 

A number of typical descent trajectories are depicted graphically. The time his- 
tories of several important trajectory parameters are presented and discussed. 


INTRODUCTION 

Evaluation of the performance capability of a specific vehicle for a lunar landing 
mission requires the selection of a basic descent mode. Possible descent modes can be 
grouped into three general categories: (1) multiple -orbit descents, (2) single -orbit 
descents, and (3) direct descents. In the multiple -orbit descent mode there is an initial 
high altitude orbit, at least one intermediate orbit at a lower altitude, and a main 
powered descent to touchdown on the lunar surface. This mode is used in the Apollo 
manned lunar landing program. If all but the final orbit are eliminated, the result is a 
single-orbit descent. The vehicle brakes directly into a low altitude orbit from which the 
main powered descent to the surface is initiated. The single -orbit mode is proposed for 
use in a lunar logistics vehicle program in reference 1. A direct descent is one in which 
the spacecraft does not enter lunar orbit prior to touchdown. The direct -descent mode 
was used in the Surveyor program and is the mode of interest in this report. 

The main advantage of a direct descent is the simplicity associated with requiring 



only one burn of the main descent propulsion system. The main disadvantage of this 
mode is that (since the perilune distance is usually less than the lunar radius) the vehicle 
approaches the moon on a collision course and will impact the lunar surface if the 
descent engine fails. For this reason the direct -descent mode is not usually considered 
for manned missions. Because of its simplicity, however, this mode is usually given 
some consideration during the preliminary planning phase of any unmanned soft lunar 
landing program. An evaluation of the vehicle’s direct -descent performance capability 
then becomes necessary. 

The characteristic velocity AV (all symbols are defined in appendix A) required to 
accomplish a direct -descent soft lunar landing at a prescribed landing site is a function 
of (1) Earth, moon, and transfer-orbit geometrical characteristics such as the transfer 
orbit vis viva energy Cg , the location of the landing site, the current Earth-Moon 
geometry and (2) vehicle propulsion characteristics such as the specific impulse I 

sp 

of the descent propulsion system and the ratio of thrust -to -vehicle weight at descent 
engine ignition (F/W). (All thrust -to -weight ratio values are based on equivalent Earth 
weight). In the preliminary planning phase of a program, the values of some of these 
variables are not well known: estimates change frequently and over relatively wide 
ranges of values, and performance analysis is time consuming. The results included 
herein are intended to facilitate the estimation of performance capability and to permit 
an evaluation of the effect of changes in many of the above variables. 

The AV requirement for direct descents is shown as a function of landing site 

location for lunar approach trajectory vis viva energy levels of 1. 78, 1. 09, and 0. 82 

kilometers squared per second squared (translunar trip times of approximately 60, 75, 

and 90 hr, respectively), descent engine I values of 300 and 440 seconds (typical of 

sp 

Earth storable and cryogenic propellants, respectively), and (F/W)^ values of 0. 12, 

0.15, 0.25, 0.50, 1.00, and 10. 0. 


ANALYSIS 

Procedure 

The mathematical problem is to determine the optimum constant thrust direct- 
descent maneuver from a lunar approach trajectory of known energy with a specified 
normal impact point (NIP) to a specified landing site on the lunar surface. The perilune 
distance of the approach trajectory and the point at which the powered descent is initiated 
are not specified and are available for optimization. The thrust direction during the 
powered descent is also unspecified and may be optimized. 

Consider the case of the specified landing site being at the NIP. Then the translunar 
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Velocity vector diagram 



Figure 1. - Schematic of lunar approach trajectory geometry. Velocity of spacecraft with 
respect to Earth at lunar sphere of influence, v ; velocity of Moon with respect to Earth, 
v em ; velocity of spacecraft with respect to Moon, v^. 


trajectory would be shaped as shown by trajectory 1 in figure 1. The spacecraft reaches 
the lunar sphere of influence with a velocity at the point P^. The velocity of the 
Moon relative to Earth v gm is subtracted vectorially from v^ to give 

Vl =7 l - ? em 

which is the spacecraft velocity relative to the moon. As shown in figure 1, v m ^ is in 
the incoming radial direction. If the constant thrust force of the vehicle is applied at the 
appropriate time and in the radial direction the vehicle will decelerate to a soft landing at 
the NIP. The perilune distance of the approach trajectory is zero in this case. 

Suppose now that the desired landing site is at the point Sg of figure 1 instead of at 
the NIP. The trajectory could again be shaped to arrive at the sphere of influence at P^, 
approach the moon along the incoming radial until descent engine ignition, and then 
employ a thrust direction profile which would allow a landing at S 2 * However, it is less 
costly, from a AV standpoint, to increase the perilune distance by reaching the sphere 
of influence at Pg and following the dashed approach path to S£ • Since the NIP is 
specified, the incoming radial of trajectory 2 is the same as that of trajectory 1. The 
asymptote of trajectory 2 is parallel to the incoming radial, and corresponds to an 
asymptote translation of magnitude d. The translunar trajectory can be designed to 
intersect the sphere of influence at Pg instead of at P^ with a velocity Vg such that 

%2 =V 2 -^em=%l 
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by making slight adjustments in the Earth injection conditions and injection time. Lunar 
approach trajectories 1 and 2 have the same energy and the same NIP but the descent AV 
requirement for landing at Sg is much lower for trajectory 2. The perilune distance of 
trajectory 2 is greater than zero. The change in perilune distance is equivalent to an 
asymptote translation of magnitude d in the direction of Sg from the NIP. The distance 
d is optimized in order to minimize the required descent AV. 

Now consider the problem of landing at the site Sg of figure 1. Site Sg is the same 
distance as Sg from the NIP but in a different (not necessarily opposite) direction. The 
minimum AV descent to Sg obviously requires a lunar approach path having an asymp- 
tote translated by the amount d from the incoming radial in the direction of Sg. The 
characteristics of the optimum descent from trajectory 3 to landing site Sg are identical 
to those of the optimum descent from trajectory 2 to site Sg (assuming a spherical, 
nonrotating moon). The minimum descent AV requirement is a function of the distance 
but not of the direction of the specified landing site from the NIP. The mathematical 
problem is therefore two-dimensional with all motion confined to the plane formed by the 
incoming radial and the approach trajectory asymptote. The descent AV requirements 
are shown as a function of the distance of the landing site from the NIP in this report. 

The descent paths are obtained by numerical integration of the two-dimensional 
descent trajectory equations of motion. For convenience, the integration starts on the 
lunar surface and terminates on an approach trajectory having the specified energy and 
NIP location. Consequently, the initial integration conditions correspond to trajectory 
conditions at touchdown, and the final integration conditions are those of the approach 
trajectory. The optimum instantaneous thrust direction xf/ and its time rate of change 
)// are determined by Lagrange multiplier calculus of variations techniques as in ref- 
erence 3. The required optimum solution and minimum AV requirement are obtained by 
solving the resulting two point boundary value problem by finite -difference Newton - 
Raphson iteration techniques. The mathematical details of the boundary value problem 
are presented in appendix B. 

The location of the NIP itself is a function of the Earth-Moon geometry, the energy 
of the translunar trajectory Cg e and the relative inclination Ai between the Moon's 
orbit plane and the trajectory plane. It is assumed that the reader will have independent 
knowledge of the location of the NIP before attempting to use the results of this report 
for performance evaluations. A procedure for estimating the NIP location follows. 


Normal Impact Point 

The AV required for a direct -descent lunar landing depends on the angular distance 
e from the NIP to the landing site and is shown as a function of e in this report. If the 
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locations of both the NIP and the desired landing site are known, e can be calculated as 


cos e = sin P n sin P f + cos P n cos P f cos Act 

where P n and Pf are the latitudes of the NIP and the landing site, respectively, and 
Act is the longitudinal difference between the NIP and the landing site. For convenience, 
the value of e is shown for all combinations of p f and Act in figure 2 for P n = 0° and 
± 10 °. 

If the NIP location is not known explicitly, it can be estimated by assuming that (1) 
the Moon moves about the Earth in a circular orbit, (2) the lunar equator is in the orbit 
plane of the Moon, (3) the Earth-Moon line intersects the lunar surface at a lunar latitude 
P and longitude ct of zero, and (4) the mass of the Moon does not have a first-order 
effect on the NIP location. These are minor approximations and provide adequate ac- 
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lb) Normal impact point latitude, ±10°. 


Figure 2. - Angular distance of landing site from 
normal impact point for all combinations of land- 
ing site latitude and longitude. 
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curacy for the type of planning estimates for which these data are intended. The NIP 
location can then be determined analytically. Details of the calculation are included as 
appendix C. The longitude and latitude of the NIP are shown in figures 3(a) and (b), 
respectively, as a function of Cg e and Ai for several values of r em . The vis-viva 
energy of the lunar approach trajectory (Cg m = v^) is given in figure 3(c) as a function 
of Co Q for the same r values. 

d) 6 Cill 


Trip Time 

Analysis of the boundary value problem and conditions shows that, for a vehicle with 
specified I and fixed (F/W)., the descent AV requirement depends only on the landing 
site location and the energy of the lunar approach trajectory. As discussed previously, 
the landing site location is specified with respect to the NIP. The energy Cg m can be 
related to the translunar trip time. The actual trip time depends on the specific tra- 
jectory, Earth-Moon geometry, etc. A good approximation of the actual trip time is 
obtained by simply calculating the two-body trip time t which is the time from Earth 
perigee passage to a radius equal to r . The value of t c is shown as a function of 

C„ for representative values of r in figure 3(d). Two-body trip time values of 60, 
o.6 em 




Ai 



(b) NIP latitude. 


Figure 3. - Normal impact point longitude and latitude, selenocentric energy and two-body 
trip time as function of transfer orbit energy for several values of Earth -Moon distance. 
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75, and 90 hours and a lunar distance of 60 Earth radii are considered representative of 

the range of interest and correspond to the values of C« used as parameters in the 

o , m 

RESULTS section. 


RESULTS 

Descent Characteristic Velocity Requirements 

The AV requirement for direct -descent lunar landings for C, „ = 1. 78 kilometers 

o , m 

squared per second squared (60 -hr translunar trip time with r em = 60 Earth radii) 

and a descent engine I of 440 seconds is shown as a function of e in figure 4(a). 

sp 

Results are included for (F/W)^ values of 0. 12, 0.15, 0.25, 0.50, 1.0, and 10.0. The 
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0 

Landing site location, deg from NIP 


60 


120 


180 


(a) Vis-viva energy (C 3 m ), 1.78 kilometers squared 
per second squared; specific impulse I sp , 

440 seconds. 


(b) Vis-viva energy C 3 m , 1.78 kilometers squared 
per second squared; Specific impulse I sp , 

300 seconds. 


Figure 4. - Descent characteristic velocity requirement as function of landing site location for various {F/W) } values. 


AV requirement at a given landing site increases rapidly with decreasing (F/W)., es- 
pecially for (F/W) i values below about 0. 5. This is a result of the increased gravity 
losses corresponding to the longer bum times for the low (F/W). values. 

For a fixed (F/W)j value the maximum AV requirement is encountered when 
landing at the NIP (e = 0). For this case the vehicle descends vertically along the incom- 
ing radial all the way from descent engine ignition to touchdown on the lunar surface. 

The thrust vector is therefore directly opposite the lunar gravity vector for the duration 
of the powered descent phase, and gravity losses are maximized. As e increases, the 
gravity losses and total AV requirement decrease. This behavior will be mentioned 
again later when the trajectory profiles are discussed. 

The AV requirement for the same trip time (60 hr) but for a descent engine I 
of 300 seconds is shown in figure 4(b) as a function of e for the same (F/W)j values as 
in figure 4(a). The characteristics of the curves in figure 4(b) are directly comparable 
to those of figure 4(a) but the actual AV requirement is lower for corresponding cases. 
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0 60 120 180 0 60 120 180 


Landing site location, « , deg from N 

(c) Vis-viva energy C 3 1.09 kilometers squared (d) Vis-viva energy C 3 m , 1.09 kilometers squared 

per second squared; Specific impulse I S p, per second squared; Specific impulse I S p, 300 sec- 

440 seconds. ond- 

Figure 4. - Continued. 

This difference is greater at the lower (F/W)^ values. The reason is that the I g p = 300 
case requires a higher propellant consumption rate than the I = 440 case for the same 

sp 

thrust and weight at descent engine ignition. As a result the vehicle weight decreases 
more rapidly and the average (F/W) during the descent is higher. At the low (F/W)^ 
values where the gravity losses are high, the increase in average (F/W) is important 
and results in reduced gravity losses. For the higher (F/W)j values gravity losses are 
small to start with and the difference in AV requirement between the 300 and 440 I 

bp 

cases is also small. 

The required descent AV for C, =1.09 square kilometers per second squared 
(75-hr two-body trip time at r gm = 60 Earth radii) and for Cg m = 0. 82 square kilo- 
meter per second squared (90 -hr two-body trip time at r gm = 60 Earth radii) is shown 
in figures 4(c) and 4(e) for I _ = 440 seconds and in figures 4(d) and 4(f) for I = 

Sp op 

300 seconds. The descent AV requirement for a specified e and a fixed (F/W). value 
is higher for the shorter trip times since these correspond to higher Cg m values. 
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Landing site location, e, deg from NIP 

(e) Trip time, 90 hours; specific impulse, 440 sec- (f) Vis-viva energy C 3 m , 0.82 kilometers squared 

onds; vis-viva energy C 3 0.82 kilometer per second squared; specific impulse I sp , 300 sec - 

squared per second squarfed; specific impulse onds. 

I S p, 440 seconds. 

Figure 4. - Concluded. 


The vehicle incidence angle at touchdown and the minimum radius prior to touchdown 
were not restricted in this study. As will be shown later the trajectories become more 
and more shallow as e increases. As e approaches 180°, the unconstrained variational 
solution produces trajectories which dip below the lunar surface. In order to obtain 
physically possible trajectories in this region it would be necessary to reformulate the 
variational equations to include a minimum radius constraint prior to touchdown. This 
was not done for this study. The AV requirement for those trajectories which include a 
subsurface dip is shown by the dashed line segments in figure 4. The dashed lines rep- 
resent a minimum bound since the constrained solution would require higher AV values. 

An actual direct descent to any landing site will require a higher AV than is in- 
dicated by the data of figure 4. No provision is included here for possible hover, trans- 
lation, and terminal descent maneuvers; no flight performance reserve has been in- 
corporated, and no mission constraints or contingencies have been included. These items 


10 




depend on the particular vehicle and mission being investigated and can be added to the 
AV requirement listed here as required. 


Trajectory Characteristics 


Descent trajectory profiles for e values between 0 and 170° are shown in polar 

coordinates in figure 5 for C, = 1.09 kilometers squared per second squared (t = 

o j m c 



75 hr at r em = 60 Earth radii), I = 440 seconds and (F/W)^ = 0. 15. The e = 0 tra- 
jectory is a vertical descent from descent engine ignition to touchdown at the NIP as 
previously discussed. As e increases the descent engine ignites at lower altitudes, 
the trajectories become more shallow, the thrust vector is more horizontal, gravity 
losses are thereby reduced, and the overall descent AV requirement decreases as was 
previously discussed. As indicated in figure 4(c) and shown in figure 5 the trajectories 
dip below the lunar surface for e values greater than about 145°. 

The effect of (F/W)j is investigated for the e = 60, I = 440, Cg m = 1.09 case in 
figure 6. The time histories of altitude, velocity, flight path angle, and thrust angle 
during the descent are shown in figures 6(a) to 6(d), respectively, for various (F/W)j 
values. Figure 6(e) shows the relation between horizontal and vertical velocity for the 
same values of (F/W)^. At the highest (F/W)j value of 10.0, the altitude and velocity 
both approach zero in nearly linear fashion and at nearly constant flight path angle and 
thrust direction values. The plot of vertical velocity as a function of horizontal velocity 
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Figure 6 . - Time history of descent trajectory parameters. Central angle between landing site and normal impact point U), 60°; specific impulse 
I n , 440 seconds; lunar approach trajectory vis-viva energy C 3 m , 1.09 kilometers squared per second squared. 



Ignition thrust 
to weight 
ratio, 



(e) Horizontal and vertical velocity component relation. 
Figure 6. -Concluded. 


is also nearly linear. For the lower (F/W)j values (longer powered descent burn times) 
the altitude decreases rapidly at first, since the vertical component of velocity is larger 
than the horizontal component, then more slowly as the time of touchdown approaches. 

As the time to touchdown is reduced the instantaneous thrust -to -weight ratio increases 
and the velocity decreases more and more rapidly. The thrust angle and flight path 
angle have their maximum value at the time of descent engine ignition, decrease to a 
minimum as the time to touchdown decreases and then increase again near landing. The 
flight path angle at time zero is not defined since the velocity is zero at that point. The 
thrust angle p at zero time is a good representation of the vehicle's attitude at touch- 
down if a tailored final descent phase is not employed. 

The effects of variations in Cg m and I on the descent profiles and parameter 
histories are relatively minor for the ranges considered in this study. 


CONCLUSIONS 

All points on the lunar surface are accessible with direct descent landing maneuvers. 
The velocity increment (AV) requirement is a strong function of the landing site location 
relative to the normal impact point and of the thrust -to -weight ratio at descent engine 
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ignition. The translunar trip time affects the direct descent AV requirement in two 
ways. First, the lunar approach velocity varies with trip time and this has a direct 
effect on the AV requirement. Secondly, the translunar trip time affects the normal 
impact point location which changes the angular distance between the normal impact 
point and the landing site and hence the AV requirement. Other factors such as the 
Earth-Moon distance and the descent engine specific impulse have a smaller effect on the 
required AV. 

The results are useful in estimating the AV required to land via the direct descent 
mode at a specified point on the lunar surface. The results can also be applied to deter- 
mine the effect on performance of changes in such things as launch vehicle capability, 
translunar injection energy or trip time, descent engine specific impulse and thrust, 
normal impact point location and desired landing site location. 

These results represent the mathematically optimum direct descent AV require- 
ments for various combinations of landing site location and lunar approach energy. The 
actual AV requirement for a specific mission must include provisions for mission con- 
straints, performance reserves and final descent and touchdown maneuvers in addition to 
the requirement shown in this report. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 5, 1970, 

731-11. 



APPENDIX A 
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SYMBOLS 


A quantity defined by eq. (Bll), m -1 

C constant of integration, dimen- 
sionless 

2 2 

C3 vis-viva energy, km /sec 

d magnitude of asymptote trans- 

lation. See fig. 1. 

n o 

E energy, m /sec 

e orbit eccentricity, dimensionless 

(F/W) ratio of thrust to vehicle equiva- 
lent Earth weight, dimensionless 

g quantity to be minimized in 

eq. (Bl), sec 

9 

h vehicle angular momentum, m /sec 

I specific impulse, sec 

Ai relative inclination between Moon's 

orbit plane and translunar tra- 
jectory plane, deg 

m mass, kg 

NIP normal impact point 

n unit normal to Moon's orbit plane 

p semilatus rectum, m 

r radius , m 

t time, sec 

AV characteristic velocity increment, 

m/sec 

v velocity, m/sec 

y flight path angle, angle between 

velocity vector and local hori- 
zontal in trajectory plane, deg 


6x allowable variation in any variable 
x 

e lunar central angle from NIP to 

landing site, deg 

ri true anomaly at descent engine 
ignition, defined in fig. 7, deg 

.6 true anomaly of incoming radial, 

defined in fig. 7, deg 

X. Lagrangian multiplier (j = 1, 2, 3, 
4, 5) 

q 9 

\i gravitational constant, m /sec 

P latitude, deg 

a longitude, deg 

Act longitude difference between land- 

ing site and NIP, deg 

0 inertial travel angle, defined in 

fig. 7, deg 

ip angle between thrust vector and 

the negative local horizontal, deg 

u> angular velocity, rad/sec 

Subscripts: 

c Earth perigee to lunar encounter 

e Earth 

•em Earth to Moon or Moon with 

respect to Earth 

f final or lunar touchdown point 

1 initial or descent engine ignition 

point 

m moon 

15 


n normal impact point 

p Earth pericenter 

°° hyperbolic 

Superscripts: 

— vector 


A unit vector 

time derivative 


* 
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APPENDIX B 


BOUNDARY VALUE PROBLEM 

As indicated in the ANALYSIS, the trajectories were integrated starting on the lunar 
surface and ending in an orbit corresponding to the lunar approach trajectory. In order 
to completely specify a vehicle's position in a two-dimensional orbit such as this, four 
independent orbital parameters must be specified. For the case at hand, only two nec- 
essary conditions are determined by the problem. The energy of the vehicle at the end 
of the integration must be equal to the energy of the lunar approach trajectory, and the 
orientation of the orbit must be such that the vehicle landing point is properly located 
relative to the NIP. The variational transversality equation can be used to find two 
auxiliary conditions which, when satisfied, insure that the unspecified orbit parameters 
assume optimum values. 

The general form of the two-dimensional transversality equation as given in ref- 
erence 3 and applied to the problem at hand is 

jc6t + X^5r + rAgSw + A^Sr + \^>(p + Ag&mJj + 6g = 0 (Bl) 

where points i and f represent the powered descent initial and final points, respec- 
tively, and g represents that quantity which is to be minimized (which, in this case, is 
the powered descent burn time). If the zero time reference is at point f (tj = 0), then L 
will be negative and 

g = t f - t. = - tj, 6g = - §t. (B2) 

The operator 6 represents the allowable variation in the variable on which it operates. 
The remaining variables in equation (Bl) are defined in appendix A and several of them 
are also indicated in figure 7. For the particular problem at hand (Bl) can be simplified 
by invoking the following considerations. 

(1) At point i the mass is specified; therefore, 

6 mi = 0 (B3) 

(2) At point f r = r m and aj = r = t = <p = 0; therefore, 

6r f = 6co f = 6f f = 6t f = 6«p f = 0 (B4) 
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(3) The mass at point f is given by 

= m- + m(t^ - t^) (B5) 

where m is the rate of change of vehicle mass during the powered descent and is 
assumed constant. Therefore, 

6m^ = + m(6tj - 6t.) = -m6tj (B6) 

Substitution of equations (B2) to (B6) into (Bl) and simplification reduce the transversality 
equation to the form: 

(C + mXg - 1)6T + Xj^Sr. + r^Xg^w. + X 3 6 i\ + X^S^ = 0 (Bla) 

All of the terms in equation (Bla) are evaluated at point i, and the subscripts will now be 
omitted for convenience. 

The two conditions determined by the problem can be used to express 5r and 6 < p in 
terms of 5r and 6c o. The first condition corresponds to the specification of the energy * 

at point i. Since energy is specified, 6 e = 0 where 

E = - (r 2 + to 2 r 2 ) - M. 

2 r 
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The variation is 


6E = r6r + cor 2 6a) + cu 2 r6r + JL 5 r = 0 

2 


which can be solved for 6r as 


6r = 



6co 



6r 


(B7) 


The second condition which concerns the orientation of the lunar approach trajectory can 
be expressed with the aid of figure 7 as 


e = 9 - rj + 


6e = 60 -6 r) + 6 <p. = 0 


(B8) 


The angle 0 is the true anomaly of the incoming radial and is expressible as 


0 = cos * — 


It then follows that 


60 = -j 6e = -V e2 -- 1 - 6p 

e|/ e 2 - 1 2e 2 p 

since, from the relations 

E = - JL (l - e 2 ) and 6 e = 0 

2p 


we can express 6e in terms of 6p as 


(B9) 


6e = e . 2 — X - 6p 
2 ep 
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The angle t] of figure 7 and equation (B8) is the true anomaly at point i and is given by 


77 = cos 


■ ■)] 


This can be written as 


cos 7] = 


■it'-) 


Taking the differential of both sides yields 


- sin tj 877 = - _j £ - 
e 




substituting 


sin 7] = 




and solving for 677 and simplifying give 




-p(e 2 + 1) - r(e 2 - 1) 


2e 2 pr 


5p + — 6r 


Substituting (B9) and (BIO) into (B8) and simplifying with the aid of the relations 


«p-i(s 

4 v 


2 4 

jup = u) r 


2o)r^5co + 4co^r^6r 




letting 


_ . Ve 2 - 1 _ w_ [p(e 2 + 1) + r(e 2 - l)jl 


r pr 


(BIO) 


(Bll) 
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and solving for 6 y make it possible to write 


5cp = 



6r + 6o> 

Ct) 


(B12) 


Substituting equations (B7) and (B12) into (Bla) and collecting terms give 



(B13) 


as the final form of the transversality equation. Since 6t, 6c o, and 6r are independent 
of each other, the only way for (B13) to be satisfied is for the coefficients of 6t, 6co, 
and 6r to all be zero. The coefficient of 6t can be zeroed by properly scaling all the 
A’s and C as explained in reference 3. The auxiliary conditions which must be satisfied 
for the descent trajectory to be optimum are that the remaining coefficients in equa- 
tion (B13) be zero. After simplification, the conditions can be written as 


-wrAj + rAg + — A^ = 0 


and 



The boundary value problem consists of finding those values of Aj to A^ at touch- 
down, the vehicle weight at touchdown and the duration of the powered descent which 
result in a trajectory satisfying the required final conditions. As in reference 3 the 
A's are scaled and initial values of Aj, A 2 , and Ag are calculated from ij/ and as 
an aid in obtaining reasonable initial guesses. Actual boundary conditions to be deter- 
mined are therefore i//, \p, A^, and vehicle weight at touchdown and burn time. The 
corresponding end conditions which must be satisfied are the specified approach tra- 
jectory energy (C 9 ), the required thrust -to -weight ratio at descent engine ignition 

O j 111 

(F/W)^, the desired angular distance from the NIP to the landing site (e), and the two 
auxiliary variational end conditions. 
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APPENDIX C 


NORMAL IMPACT POINT DETERMINATION 


If the four assumptions listed in the ANALYSIS are satisfied then the 
can be determined analytically as follows. For specified values of C 3 e 
velocity v and flight path angle y of the spacecraft at a radius equal to 
calculated as (see ref. 5): 


NIP 

location 

and 

r . the 


em’ 

r 

em 

can be 


em, 


1/2 


( 


r 

cos y = — 2- 


em 


c^i 

em 


where p is the gravitational constant of the Earth and r p is the Earth perigee radius 
of the translunar trajectory. The velocity v along with the angles y and Ai completely 
specify the vector velocity v at lunar intercept as shown in figure 8. The relative in- 
clination Ai is measured from the Moon’s velocity vector v em to the local horizontal 
(where the local horizontal is in the trajectory plane) and can have any value from 0 to 
360°. 
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Figure 8. - Velocity vector diagram at lunar 
orbit intercept. 



The Moon’s velocity relative to Earth v gm is orthogonal to the Earth -Moon line because 
of the circular orbit assumption and has magnitude 



The spacecraft velocity relative to the Moon is 


The direction of v tc defines the incoming radial direction. The NIP is the intersection 
of the incoming radial (or of -vj and the lunar surface as shown in figure 9. 



Figure 9. - Normal impact point location. 
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